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Abstract. The hot, puffed up, post-shock region of an advective disc is the source of 
high energy photons and also the jets and outflows. We study the relativistic equations of 
motion of jets as these high energy photons interact with them. We show that the much 
discussed terminal velocity of jets depends on the comparative value of radiative energy 
density, flux and the radiative pressure. We show that electron-positron pair plasma jets 
achieves highly relativistic terminal speeds for higher disc luminosities. 



1. Introduction 

Rotating matter while falling onto black holes creates a temporary depository of 
matter called accretion disc. As matter moves closer to the black hole, at around 
a few tens of Schwarzschild radii (r g ), centrifugal force tends to be comparable to 
the inward gravitational force and the supersonic inflow is slowed down. If this 
slowing down occurs in a thin region [1], the flow is said to suffer a shock, where 
the Mach number jumps discontinuosly from supersonic to subsonic branch. This 
slowed down post-shock flow is considerably heated up and as a result, puffs up 
(to maintain hydrostatic balance along z-direction) in the form of a tori called 
CENBOL (centrifugal pressure supported boundary layer). CENBOL contains a 
copius hot electrons which inverse-Comptonize the soft photons from outer cool 
thin disc [2] , producing the hard-power law tail [3] . As a lot of heat is stored in 
the CENBOL, Chakrabarti and his co-workers (e.g. [4], [5], [6]) has shown that 
the thermal pressure along the vertical direction pushes the matter out in the 
form of jets. In this paper, we investigate the interaction of these jets with the 
high energy photons radiated by the CENBOL. 

Interaction of radiation and jets have been studied by several workers. Piran 
([7]) showed that it is difficult to accelerate outflows from thick accretion disc 
beyond Lorentz factor 7 > 1.5. Icke ([8]) in a very important paper showed that, 
if outflows are to be accelerated by radiation from 'infinite' thin discs, the terminal 
speed ('magic speed' in his parlance) achieved is just 0.451c, where c is the velocity 
of light. Fukue ([9]), showed the terminal speed for rotating flow above thin disc, is 
little less than what Icke had obtained. The conclusion is that radiative interaction 
actually limits the jet terminal speed. We improve these works by studying the 
interaction of a jet with radiation coming from the most general form of accretion 
flow, namely the advective discs (e.g. [1],[3],[10]). The first three moments of the 
radiation intensity for e.g. , the radiation energy density, radiative flux and various 
components of radiative pressure are being calculated following the treatment of 
Chattopadhyay & Chakrabarti (e.g. [11], [12]). For simplicity, we assume that the 
radiation originates only from the CENBOL. 
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In the next section we show that, the terminal speed achieved by the outflow- 
ing plasma depends on the comparative value of the various moments of radiation 
involved. We show that the radiation from the advective flows can produce rela- 
tivistic jets. In §3, we draw our conclusions. 
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Figure 1. Comparison of radiation field quantities £ (solid), T (dashed), and V (long-dashed) 
plotted with log(z) in non-dimensional units. These field quantities are calculated for disc lu- 
minosity L — OML-Edd- Radial dimension of CENBOL (i.e., shock location) is chosen to be 
Wr g . 



2. Equations of motion and radiative acceleration of jets 

We assume that the electron-positron plasma jets are non-rotating and con- 
fined along the axis of symmetry. To study only the radiation effects, we ig- 
nore the gas pressure gradient term. We confine our investigation within the 
realm of special relativity and the effects of strong gravity is taken care by 
Paczyhski-Wiita potential (e.g. [13]). The unit of length, time, mass and velocity 
are r g = 2GMb/c 2 , 2GMb/c 3 , Mb and c respectively, where G and Mb are the 
gravitational constant and mass of the central black hole. The metric considered 
is given by, dr 2 = dt 2 — dr 2 — r 2 d(j) 2 — dz 2 . We assume u r = = and also 
d/dt = d/dr = d/d(j) = 0. We define a 3-velocity v, given by v 2 = —(uiU l )/(utu t ), 
where the latin index (e.g. i,j,k) denotes space variables. The Lorentz factor is 
given by, 7 2 = 1/(1 — v 2 ). Hence in our case u z = jv, u z = — jv, u l = u t = 7. 
The equation of motion (see, Mihalas and Mihalas [14] and Kato et al. [15]) is 
given by, 

dn z 1 

"JF = ~ 2(z-l) 2 + Kes[lFZ ~ ^ EuZ ~ PZZuZ + uZ ^ FZuZ ~ pZZ u z u z )}, (1) 
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where, K es , F z , E, P zz are the Thomson scattering opacity, the z-component of 
radiative flux on the z-axis, energy density on the axis of symmetry, and the z-z 
component of pressure tensor on the axis of symmetry. The term — l/{2(z — l) 2 } 
is the gravity term with Paczyhski-Wiita potential. Using the expressions of u z 
and u l we can reduce Eq. (1) in terms of the 3- velocity v, and is given by, 



dv_ -2(1^ + ~ l 3 v£ - TvP_ + 7 3 (2^ ~ v 3 V)} 
dz 7 4 i> 



(2) 



where, J 7 , £, and V are radiative flux, energy density and z-z component of 
pressure tensor which has been multiplied by n es . To find the expression for 
terminal speed, we put the term in the square bracket of Eq. (2) equal to zero 
and get a quadratic equation for v t , 

Tv\ -{£ + V)v t + F = 0. 

It is easy to show that at z^large, V ^ T ^ £, and also that all three quantities 
vary slowly with z (Fig. 1). Thus for large z, let V = A, hence T = A + 5 and 
£ = A + i], where S-^A and r]<^A with S ^rj. Hence, from the above quadratic 
equation we get, 

„ (£ + V)-V(£ + vr-^ 

Vt 2T (3) 

Putting the radiation field quantities in terms of A in Eq. (4), we have, 

Equation (4) shows that vt can be very close to the velocity of light, thus radiation 
drag does not limit the terminal velocity to some moderately relativistic values. 
In case of gas flows above a thin disc (e.g. [8], [9]), V = 1/35 and £~2.F, putting 
those values in Eq. (4), we have vt = |(4 — y/7) (see, Icke[8]). We have thus 
proved that the terminal velocity achieved depends on the radiation field produced 
by the structure of the disc. If we now integrate Eq. (2) with a very small injection 
velocity, we can find the solution topologies for jets as shown in Fig. 2a. Though 
there is no upper limit to vt, even then increasing the luminosity, will not increase 
Vt drastically, as there is 7 4 term in the RHS of the denominator of Eq. (2) , and 
hence as v— >1, dv/dz^O. In Fig. 2b the velocity at z = 10000r s i.e. , wioooo 
is plotted with the disc luminosity L and is found that as one increases L, the 
response of wioooo is not linear. Figure 2b shows that it is difficult to get terminal 
speeds ~0.99c with radiative acceleration only, as that would require extremely 
high luminosity, which is difficult to produce without, at the same time, cooling 
the CENBOL itself. 



3. Discussion and Concluding Remarks 

In our investigation we have assumed the entire disc luminosity is coming only 
from the CENBOL region. It can be shown elsewhere that the contribution from 
the pre-CENBOL thin disc will marginally affect the results shown here. We con- 
clude that, the terminal speeds calculated previously, by using radiation fields 
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Figure 2. (a) Velocity profile of an electron-positron jet, with input parameters Zi„ = 3r s , 
Vin — 0.05c. The radiation field quantities correspond to disc luminosity 0.061/Edd- (b) The 
variation of Wioooo {i- e - v at z = 10000r 9 ) with the accretion disc luminosity L in units of LEdd- 



from very specialised accretion discs {e.g. , thin disc, thick disc, etc), show mildly 
relativistic flows, and are basically the artefacts of those disc geometries them- 
selves. If one considers more generic disc model as we have done, the radiation 
drag does not introduce a mildly relativistic upper limit for terminal speed. We 
find that achieving terminal velocity close to that of light by radiation pressure 
effects is possible. 
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